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Abstract 

Although the clean Si(001) and Ge(001) surfaces are very similar, experiments 
to date have shown that dimer-vacancy (DV) defects self-organize into vacancy 
lines (VLs) on Si(001), but not on Ge(001). In this paper, we perform empirical- 
potential calculations aimed at understanding the differences between the vacancies 
on Si(001) and Ge(001). We identify three energetic parameters which characterize 
the DVs on the two surfaces: the formation energy of single DVs, the attraction 
between two DVs in adjacent dimer rows, and the strain sensitivity of the formation 
energy of DVs and VLs. At the empirical level of treatment of the atomic interac- 
tions (Tersoff potentials), all three parameters are favorable for the self-assembly 
of DVs on the Si(001) surface, but not on Ge(001). The most significant difference 
between the defects on Si(001) and on Ge(001) concerns the formation energy of 
single DVs, which is three times larger in the latter case. By calculating the strain- 
dependent formation energies of DVs and VLs, we propose that the experimental 
observation of self-assembly of vacancies on clean Ge(001) could be achieved by 
applying compressive strains of the order of 2%. 

Keywords: Surface defects, Self-assembly, Surface energy, Semi-empirical models 
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1 Introduction 



Epitaxial growth of germanium on the Si(OOl) surface has been an active area of investi- 
gation because it represents a versatile way of fabricating quantum dots. At the same time, 
Ge/Si(001) provides an example of a system where the mismatch strain can be used to control 
the morphological features of the surface. In the early stages of growth, the surface changes its 
reconstruction from the 2x1 structure to a 2 x N pattern, which consists of an array of long, 
parallel lines of dimer vacancies. It has recently been shown that the 2 x N reconstruction can 
be an excellent template for the nanofabrication of 1-dimensional structures with unique elec- 
tronic transport behavior 0[21 El- At later stages of growth, the film breaks up into pyramidal 
islands (quantum dots); controlling the size, spatial distribution and electronic properties of 
these quantum dots represent key issues for the practical applications of the heteroepitaxial 
growth on Si(OOl). From a fundamental point of view, the current state-of-the-art in scanning 
tunnelling microscopy has helped renew the general interest in the 2xN structures on Si(OOl). 
For example, Rastelli and coworkers [I] have recently discovered that the 2 x N reconstruction 
on the Ge-covered Si(OOl) surface disappears upon Si capping, confirming in an ingenious way 
the role played by the mismatch strain in determining the stability of the vacancy lines. In 
another recent study 0, Sutter et al. have used high-resolution STM imaging in conjunction 
with a detailed statistical analysis of morphological features of the surface to establish that 
the interactions between steps and vacancy lines cause the initial surface roughening of the 
Ge/Si(001). 

Rapid progress has been made after the discovery of the 2 x N structure 00, and two 
authoritative reviews summarizing the present understanding of the growth of germanium 
films on Si(OOl) have been published Dimer vacancies (DVs) form easily on the Si(OOl) 

surface, but they only organize into vacancy lines (VLs) when there is a sufficient amount of 
Germanium deposited (close to 1 ML). The mismatch strain created by epitaxial deposition 
increases the number of dimer vacancies per unit area, and the interactions between these 
vacancies become important at sufficiently high vacancy concentrations. The interaction be- 
tween two DVs is long-range repulsive when the DVs belong to the same dimer row, but 
weakly attractive when they belong to adjacent dimer rows lllj. Such anisotropy of the 
interactions determines the alignment of the DVs in the direction perpendicular to the dimer 
rows. This physical picture for the origin of the 2 x N reconstruction has emerged from the 
experimental work of Chen et al. |lUj . The mechanism for the formation of vacancy lines is 
consistent not only with experiments done on other epitaxial systems (e.g., Bi/Si(001) 12 j. 
but also with the observations of the Si(OOl) surface with metal impurities such as Ag ^3] or 

NidUCE]. 

While the physical origin of the vacancy lines on Si(OOl) has been elucidated [TJU, it 
remains a puzzle as to why vacancy lines are not normally observed on the very similar 
Ge(OOl) surface. The fact that the Ge(OOl) surface is virtually defect free |E3 QZ| suggests 
that the formation energy of the vacancies on Ge(OOl) is substantially larger than that of the 
vacancies on Si(OOl). Quantitative data from atomic-scale calculations on the Ge system is 
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not available, although the need for such calculations was emphasized almost a decade ago in 
the work of Yang et al. 18 . The purpose of our paper is to address the differences between 
Ge(OOl) and Si(OOl) concerning the formation of dimer-vacancies and their self-assembly into 
vacancy lines. The work presented here fits well into the current efforts for explaining 
the differences between the two technologically important surfaces. 

We have calculated the formation energies of point and line defects on the (001) surfaces 
and their dependence on applied strain, using empirical potential |19| 120 j models for interac- 
tions . The results obtained offer an explanation for the differences between the vacancies on 
Si(001) and Ge(001) and are consistent with experimental findings; moreover, these results 
could foster further work based on more accurate descriptions of the atomic interactions such 
as tight-binding or density functionals. We have considered both the heteroepitaxial system 
Ge/Si(001) for different values of Ge coverage, as well as pure Si and Ge surfaces under exter- 
nal strain. While the results concerning the epitaxial Ge/Si(001) systems have been previously 
obtained in different forms |21l 1221 OH] , the variation of the energetic parameters of DVs and 
VLs on the (001) surface with applied strains has not been systematically addressed in the 
past. Such a study of the energetics of defects on (001) surfaces as a function of applied strain 
is presented here. The new results obtained in this paper are: 

(a) the formation energy of the dimer-vacancies on Ge(001) and its decrease under compres- 
sive, biaxial strains. In the case of unstrained surfaces, we have found that the formation 
energy for single DV on Ge(001) is three times larger than the corresponding value for the 
single DV on Si(001). Therefore, in order for the formation energy of a single- DV on Ge(001) 
to become negative, a large compressive strain should be applied, either biaxially or in the 
direction of the dimer rows. 

(b) the formation energy of vacancy lines on Si(001) and Ge(001) as a function of applied 
strain. Using a simple model based on the interactions of nearest-neighbor vacancies in a 
VL, we have also identified the formation energy of a single DV and the attractive interac- 
tion (binding) between two DVs in adjacent dimer-rows as the dominant contributions to the 
formation energy of a VL. 

The organization of this article is as follows. In section 2 we outline the competition 
between the formation energy of the vacancy lines and their elastic interactions as the general 
energetic mechanism that leads to the appearance of the 2 x N reconstructions. In contrast to 
previous interpretations [S], we have found that there is no long range attraction between the 
vacancy lines, as their interactions are purely repulsive. As in most other studies to date, the 
formation of vacancy lines is regarded here as a purely thermodynamic phenomenon, thus the 
kinetics concerning the diffusion of vacancies during their alignment are not considered. In 
section 3 we describe the computational details of the total energy calculations based on Tersoff 
potentials |19[ I20j . The next two sections are devoted to the formation energies and repulsion 
strengths for single-vacancy lines (VLs) on Ge-covered Si(001) (section 4) and for different 
types of VLs on pure Si(001) and Ge(001) (section 5). We investigate the formation energies 
of individual dimer vacancies on (001) surfaces and their strain dependence, identifying the 
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differences between the defects on Si and Ge surfaces. The possible relevance of these results 
for future experimental work is also discussed in section 5. Our conclusions are summarized 
in section 6. 



2 Energy of vacancy lines on 2 x N surfaces 

The 2xN reconstruction is a structural pattern on the (001) surface of Si, characterized by the 
periodic lengths 2a and Na, where N is a positive integer and a = 3.84A is the lattice constant 
of the unreconstructed Si(001) surface. This reconstruction is obtained by eliminating every 
iV-th dimer from each dimer row, such that the dimer vacancies created in this manner are 
in registry with one another and form straight lines perpendicular to the dimer rows. To 
lower the number of dangling bonds at the surface, the second-layer atoms below the dimer 
vacancy rebond such that they are fully coordinated, at the cost of introducing some surface 
stress [^|22l|2Hj. Depending on the preparation details of the surface, other types of atomic 
structures of VLs can be obtained El El El El HH ESI HHj- in this section we focus on 
the general features of the energetics of the 2 x N surface, without reference to a specific 
atomic structure. We start by partitioning the surface energy into a contribution from the 
2x1 surface and a contribution from the vacancy lines: 

A 

72x7V = 72x1 + (1) 

where 72x7V and 72x1 denote the surface energies of 2 x N and 2x1, respectively, A is the 
energy per unit length of the vacancy lines, and 1/L X is the density of lines in a periodic 
array of VLs with a spacing of L x = Na . Eq. Q defines the energy per unit length A of the 
vacancy lines as an excess energy with respect to the defect-free 2x1 surface. We note that 
the energy of the vacancy lines A includes the formation energy of the lines, as well as their 
elastic interactions. Guided by the early work of Marchenko and Parshin [22j, we expect the 
interactions to be inversely proportional with the square of the separation between VLs, so 
that the line energy A can be written as 

where A is the formation energy of a VL and G is the strength of the repulsion between lines. 
The numerical factor 7r 2 /6 arises since we treat a periodic array of line defects |28j . rather 
than two isolated VLs. Using Eq. (J2J) in (^Q) we can write the surface energy of the 2 x N 
reconstruction as a function of N: 

A 7T 2 G 

where we have used the fact that L x = Na. 
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Figure 1: A typical plot of the surface energy difference 72x/v — 72xi as a function of N (solid 
diamonds). The minimum of the surface energy 72x/v is determined by the competition be- 
tween the negative formation energy (triangles) and a repulsive (positive) interaction between 
the vacancy lines (circles). 

Analyzing Eq. ©, we note that: 

(a) if A is positive then 72x7V > 72x1, so the 2x1 structure is energetically favorable. 

(b) for any A < 0, Eq. © has only one minimum located at: 



For the case of negative formation energy A, a typical plot of the surface energy difference 
(72x/v — 72x1) as a function of N is shown in Fig.^ The condition A < is sufficient to ensure 
the stability of the 2xN reconstruction (over the 2x1 structure), because for N = N* we have 
72x7V — 72xi = (2/3)A < 0. This analysis of Eq. © shows that, irrespective of the atomic-scale 
features of the VLs, the 2 x N reconstruction is always stable when the formation energy of 
a vacancy line is negative. The optimal value of N is given by the competition between the 
formation energy A and the strength G of the repulsive interactions, as shown in Eq. Ijljl. The 
magnitudes of A and G are determined by atomic-scale details, and incorporate the effects of 
the VL structure [2D 122] , wetting layer composition ESI EHj , as well as the anisotropy of 
surface stress and relaxation |22l 12*3] . 

The following two sections are devoted to the computation of the line formation energies 
and repulsion strengths for various VL structures, and to the investigation of the behavior of 
these quantities as a function of germanium coverage and applied strain. Before we proceed 




(4) 
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with the numerical results, we draw a connection with previous work on the energetics of 
the 2 x N surface. The above theory of the energetics of the vacancy lines is similar to 
Refs. |3U1 131j . since it explicitly considers the competition between the (negative) formation 
energy of a VL and the repulsive (positive) interaction between the VLs as the physical origin 
of the 2xN structure. In Ref. |30| however, the emphasis falls on the surface stress anisotropy, 
and the treatment of the two-domain (i.e. 2 x N and N x 2) surfaces obliterates the simplicity 
of the VL energetics outlined above. 

A word of caution is due concerning the mechanisms of formation of the 2 x N recon- 
struction proposed in Ref. (2H1 and in the review article [H] ■ In [21]] j the conclusions based on 
the energetic competition between "trench formation" and "strain relaxation" are fortuitous, 
since the TV-dependence of each of these two energetic contributions is incorrect . The sign 
of the "trench energy" and the magnitudes of the two energetic contributions are also incor- 
rect, though these aspects are probably outside the scope of the simple Frenkel-Kontorova 
model used in "2H] . Ref. [S] asserts that "in the dimer-row direction the VL-VL interaction 
is short-range repulsive and long-range attractive", which seems plausible given the behavior 
of the 72x7V in the limit of large N. However, we note that what appears to be an attractive 
long-range interaction between VLs is in fact the formation energy of one VL per unit cell of 
the 2 x N surface. As shown in the sections to follow, the VL-VL interaction takes a purely 
repulsive form once the vacancy line energy is extracted from the surface energy. 

3 Computational Details 

3.1 Structure of the vacancy lines 

In the terminology of Ref. [22], the vacancy lines considered here are made of 1-DV, 2- 
DV and 1+2-DVs, as shown in Fig. [21 The choice of these structures is motivated by the 
experimental observations of the atomic-scale details of vacancy lines on different surfaces: 
the 1-DV structure is present on Ge-covered Si(OOl) [HI O [221 EH] , while lines with the 2- 
DV and 1+2-DV configurations have been observed on Si(OOl) with Ni [°"IJ El and Ag 
impurities, as well as on radiation-quenched Si(OOl) [2*11 [23j . 

3.2 Total energy calculations 

Accurate calculations of the long-range behavior of the VL-VL interactions are very demand- 
ing, requiring slabs that are both long (N > 10) and thick (at least 10 layers). Earlier density- 
functional and tight-binding calculations were restricted to thin slabs with 4-6 moving layers, 
due to limitations in computational power |33[ I34 | l35j . At present, density- functional and 
tight-binding calculations have become possible |36M37| . though they remain computationally 
expensive for the purpose of determining the long-range interactions between line defects. 

1 We believe that the energetic competition that determines the 2 x N reconstruction is described by Eq. J5J). 
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Figure 2: Atomic structures of the vacancy lines, with the nomenclature given in Ref. |32| . 
The 1+2-DV complex is made of a 2-DV complex and single DV, separated by one dimer row. 

Such calculations are still prohibitive for point defects on the surface, since these defects re- 
quire simulation cells that are large in all three directions. As more studies would probably 
become available for various types of vacancy lines in the near future, it is worth pointing 
out that a good test for their total-energy convergence is to check that the line energies A 
can be fitted well with the inverse-square power law given by Eq. (j2j). Such a test would also 
indicate whether or not the chosen thickness of the computational slab is sufficient. This type 
of analysis is routinely performed for steps on crystal surfaces (see |381 1391 14U| and references 
therein) in order to determine the step formation energy and step-step interactions, but it has 
almost never been applied to the case of vacancy lines. 

Since we will investigate vacancy lines as well as individual vacancies (point defects) 
on the surface, we employ empirical potentials to model the atomic interactions. We use 
the Tersoff potential for Si and Ge ^HEO], since it satisfactorily describes the structure and 
energetics of steps on (001) surfaces [23 SH, as weu as the surface energies of the 2 x N 
reconstructions [2^]. Despite known limitations, the empirical potential described in |20j 
represents a widely accepted standard for Si-Ge systems with large numbers of atoms. Our 
aim is restricted to a meaningful comparison (rather than an accurate determination) of the 
properties of dimer vacancies on Si(001) and Ge(001). For this purpose, the use of the Tersoff 
potential is appropriate, and even desirable as a starting point. 

For simulation cells with a thickness of 200A and values of N in the range 4 < N < 36, 
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we have performed conjugate- gradient structural optimizations. From the minimized total 
energy E of the atoms in the simulation cell we compute the surface energy 72 xN- When 
there is only one atomic species in the simulation cell (e.g., Si), the surface energy is given by: 

E = ne S i + 72xNL x L y , (5) 

where esi is the bulk cohesion energy per Si atom and n denotes the number of atoms in the 
simulation cell that are allowed to relax. For the case of two atomic species in the simulation 
cell, the surface energy is given by: 

E = n Ge fi Ge + nsi^si + 12xnL x L v , (6) 

where na e , H-Ge {nsi-> (J-Si) are the number of atoms and chemical potential of Ge (Si). Since we 
are considering deposition on Si(OOl), the chemical potential of Si is set to the bulk cohesion 
energy per atom esi- For low Ge coverages (1-3ML), the chemical potential of Ge (nGe) 
depends on the deposition conditions and is not known, although its value is usually set 
to the bulk cohesion energy of (unstrained) Ge, ec e [HI H31 EH ESS • Exceptions are found 
in the work of Oviedo et al. [HE] and Li et al. (SZj, who carefully define a reference state 
from which the surface energies are measured; the procedure in |36l 13 7j is equivalent to 
computing the chemical potential of Ge based on a wetting layer composition assumed to 
be thermodynamically stable. From the results of Refs. |36l I37j we infer that for a pure Ge 
wetting layer (no intermixing), setting [iQ e = ec e is a poor approximation for monolayer 
and sub-monolayer coverages, but it becomes reasonable for 2 ML and 3 ML of Ge coverage. 
Though we use nee = &Ge in the present work, we will also address the effects of changing 
the chemical potential of Ge. 

4 Vacancy lines on the Ge-covered Si(OOl) surface 

In this section, we calculate the surface energy of the 2 x N surface with 0-3 ML Ge coverage, 
then compute the line formation energy and the VL-VL interactions. Our results for the 
surface energies are very close to the ones originally obtained by Liu and Lagally [22]; the 
deviations are due to the slightly different versions of the Tersoff potentials 2 . What differen- 
tiates the results of this section from those of Ref. [23] is that we pursue the determination of 
the line formation energy and VL-VL interaction strength, whereas Liu and Lagally address 
the interplay between surface stress and stoichiometry |23| . 

Using Eq. (^Q), we have computed the line energy A and plotted it as a function of the 
inverse square of the line spacing (l/L^) in Fig. [21 The line energy A is a linear function 
of l/L^ for each integer value of the Ge coverage between OML and 3ML (Fig. EJ). This 

2 The parameter A3 of the Tersoff potential given in |19| is set to zero in |2()|. thus modifying the original T3 
parametrization for pure Si. Here, we have found it convenient to use a form of the interactions that reduces 
to |19| in the limit where only one atomic species is present in the simulation cell. 
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Figure 3: Vacancy line energy A plotted as function of inverse square of the line spacing 
(1/L|) for different Ge coverages: ML (diamonds), 1 ML (circles), 2ML (triangles) and 3 
ML (squares). The inset represents the energy of the 2 x N surface as a function of N |23| . 
computed with the Tersoff potential T3 [T§1 [2*0] . 

behavior confirms the expectation (Eq. (J2J) ) that the interactions between the vacancy lines 
are dipolar in nature [21] • The parameters of the linear fits shown in Fig. |H1 give the line 
formation energy A (intercept) and the line- line repulsion strength G (slope), and are listed 
in Tabled The table also includes the fitting parameters for two cases of intermixed wetting 
layer compositions. We have verified that the the dipolar character of the line-line interactions 
is not changed by intermixing in the wetting layer. These intermixing cases are considered 
here only for the purpose of illustration, as the actual structure of the wetting layer at given 
Ge coverage is under debate [H E3 E2 E2] . The data in Table □ shows that for both of the 
wetting layer compositions considered (pure Ge or intermixed), the formation energy decreases 
strongly with increasing Ge coverage, which is consistent with previous findings (refer to jHj). 
The repulsion strength G is very sensitive to the WL composition, as the elastic interactions 
are largely transmitted through the first few layers: different intermixed compositions affect 
the elastic properties of the film and thus the strength of the repulsive interactions. 

The trends shown in Tableware well understood |Hj and physically meaningful, although 
we emphasize that the exact values of the formation energy A are influenced by the choice 
of the chemical potential for the deposited Ge atoms. By taking finite differences in Eq. l|S|). 
we note that a change 5fj,Ge in the chemical potential of Ge introduces a change d^xN in the 
surface energy of the 2 x N structure given by: 

= n Ge (5fi Ge ) + (2N a 2 ) (S l2xN ), (7) 
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WL 


A(meV/A) 


G (eVA) 


N* 


ML 


-20.44 


15.562 




1 ML 


-27.95 


15.603 


14 


2 ML 


-129.92 


13.526 


6 


3 ML 


-134.61 


13.818 


6 


1.5 ML (i) 


-11.03 


7.968 


15 


2 ML (i) 


-12.20 


8.541 


15 



Table 1: Formation energy A of the vacancy lines and their repulsion strength G calculated 
with the Tersoff potential [20] for the Ge-covered Si(001) surface. The first four rows cor- 
respond to a wetting layer (WL) of pure Ge 0-3 ML thick, and the last lines represent two 
intermixed (i) structures with the composition given in Ref. |23| . The last column shows the 
values of A^* obtained from Eq. @. 



where we have used the values for the periodic lengths, L x = Na and L y = 2a. Similarly, for 
a 2 x 1 surface with the same area, thickness and WL composition, we have 

= (2 + n Ge ){5fi Ge ) + {2Na 2 ){5 l2 xi), (8) 

where the number of Ge atoms was increased by 2, as no dimers are missing from the 2x1 
structure. Subtracting Eq. (|HJ) from Eq. (jJJ) and using Eq. (J2J) we obtain 

4 , Gt = A ^(_ + __), (9 ) 

where 5 A and 5G are the changes in the line formation energy and repulsion strength generated 
by the change in chemical potential 5fj.ce- From Eq. @ in the limit of large N, we find that 
8A = 8/iGe/o. and 5G = 0. Therefore, modifications of the Ge chemical potential amount 
to shifting the line formation energy while leaving the VL-VL interactions unchanged -a 
conclusion which is independent of the model used to describe the atomic interactions. 

In summary, we have computed the line formation energies and the repulsion strengths 
for certain Ge coverages and WL compositions. As shown in TableQ the line formation energy 
A is sensitive to the WL thickness and composition, while the VL-VL interactions are mostly 
sensitive to composition. To gain further insight into how strain affects the interactions and 
self-organization of vacancy lines, we will now consider 2 x A^-reconstructed surfaces of pure 
species (either Si or Ge) under external applied strain. 



5 Defects on Si(OOl) and Ge(OOl) under applied strain 

The study of vacancy lines and dimer vacancies on uniformly strained Si(001) and Ge(001) 
surfaces has two technical advantages: (a) the ambiguity regarding the value of chemical 
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VL 

structure 


Go 
eVA 


A 

meV/A 


uniax (x) 


A'(meV/A) 
uniax (y) 


equibiax. 




1-DV 


15.56 


-20.56 


4868.5 


-412.55 


4468.5 


Si 


2-DV 


5.01 


51.57 


2879.3 


-822.08 


2067.7 




1+2-DV 


33.33 


30.56 


6530.6 


-1228.5 


5354.2 


Ge 


1- DV 

2- DV 
1+2-DV 


11.12 
4.00 
24.36 


32.48 
69.613 
104.64 


4112.0 
2543.1 
5617.1 


-559.72 
-846.11 
-1396.5 


3559.3 
1709.6 
4230.8 



Table 2: Repulsion strength Go and line formation energy Ao at zero strain for the vacancy 
line structures shown in Fig. |2J The last three columns show the strain sensitivity A' of the 
line formation energy for different types of in-plane deformation. The results in this table are 
obtained using the Tersoff potentials ^3 [2U] . 

potential of Ge is removed, and (b) the uncertainties related to intermixing in the wetting 
layer obviously do not arise when studying systems with only one atomic species. Therefore, 
the study of defects on clean surfaces under external strain allows for a more direct comparison 
between Si(001) and Ge(001). 

5.1 Vacancy lines on strained surfaces 

The surface energy of the 2xN structure under external strain can be computed using formula 
J5J) in which the bulk cohesion energy is also a function of the applied strain. The in-plane 
deformations are applied along the dimer rows (e x = e, e y = 0), perpendicular to the dimer 
rows (e x = 0, e y = e), or equibiaxially (e x = e y = e). For each type of deformation we vary 
the strain in the range —2% < e < 2%, by appropriately scaling the atomic coordinates. In 
this range of deformation, we find that the line formation energy A is a linear function of 
strain: 

A{e x , e y ) = A + A' a e, (10) 

where Ao is the line formation energy in the absence of strain, and the subscript a denotes the 
uniaxial strain along different directions (a = x, or a = y) or the equibiaxial strain (a = xy). 
A consequence of the linear relation described in Eq. (jlOj) is that the biaxial strain sensitivity 
A' xy is the sum of the two sensitivities A' x and A' y that correspond to the uniaxial deformations: 

A' xy = A' X + A' y . (11) 

For the VL structures shown in Fig. |2 we have computed all the strain sensitivities 
A' a (a = x, y, xy) independently, and found that formula (fTT|) holds to within 1-2% of 
A' xy . These strain sensitivities are listed in Table |2J along with the formation energies Ao 
and repulsion strengths Gq computed in the absence of strain. The strain-dependence of the 
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line formation energy A is strongly anisotropic. As shown in Table |2 the strain sensitivity 
A' x along the dimer row is positive, while A' y is negative and several times smaller (up to one 
order of magnitude for 1-DV vacancy lines) in magnitude than A^. This indicates that for 
all the vacancy structures shown in Fig. |^1 the rebonding in the second layer plays a crucial 
role in determining the strain dependence of the line formation energy. In the absence of 
strain, the bonds that "bridge" the second-layer atoms in the vacancy lines are stretched 
compared to their bulk values, and are especially sensitive to strains applied parallel to them. 
Therefore, compressive strains along the dimer rows will decrease the amount of stretch in 
those bonds, lowering the formation energy of the vacancy lines. This correlation between 
atomic bonding and formation energy is consistent with the ones given previously for the case 
of 1-DV lines (see, for example |8l l21[ l22 l l23|). On the other hand, an uniaxial strain applied 
perpendicular to the dimer rows only indirectly affects the formation energy of the vacancy 
lines via a Poisson-type effect: a stretch along the y-direction creates a small compressive 
stress in the x-direction, which in turn lowers the line formation energy. 

While the formation energy A of the vacancy lines varies significantly with applied strain, 
the VL-VL repulsion strength G shows a weak strain dependence, changing by less than ten 
percent for the values of strain considered here. The reason for the weak strain-dependence 
of the VL interactions is that the relaxations that occur after applying external strains to 
an already relaxed 2 x N surface are rather small, since no bonds are broken or formed. 
Consequently, the VL-VL interactions are not substantially affected by external strains as 
long as these strains remain in the small-deformation regime given by Eq. (|10l) . 

By comparing the Si and Ge data (Table [2j) for each type of VL structure, we find 
several systematic trends: (a) the repulsion strengths G are stronger for Si than for Ge, (b) 
the formation energies Ao are smaller for Si than for Ge and (c) the strain sensitivities along 
dimer row (A x ) and biaxial (A xy ) are larger for Si than for Ge. Most notably, the formation 
energy of a single-DV line on Si(OOl) is negative, while for the case of Ge(OOl) it is positive. 
While there is a possibility that this sign difference between Si(OOl) and Ge(OOl) might be due 
to artifacts of the empirical potential, it is consistent with the observation of self assembly of 
DVs on Si(OOl) 26, and not on Ge(OOl) [TBI EEB1 ESI EH • Some of these trends are related to 
the formation and binding of individual dimer vacancies, as will be shown in the next section. 

The occurrence of the vacancy lines is a consequence of the strongly anisotropic inter- 
actions between the individual dimer vacancies, coupled with the preference of the DVs to 
diffuse along the dimer rows. As discovered by Chen et al. |10| . two DVs that belong to 
the same dimer row repel one another, while DVs in adjacent dimer rows experience a weak 
attraction at short distances. This finding was confirmed by other independent experiments 
[141 113j , as well as by theoretical studies . Thus the physics of vacancy lines is intrinsically 
related to the formation, interactions and diffusion of individual dimer- vacancies. In what 
follows, we show how the strain dependence of the formation energy of the VLs is linked to 
the behavior of individual dimer vacancies under strain. 
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5.2 Dimer vacancies on Si(OOl) and Ge(OOl) 

In this section, we compute the formation energy of three types of DVs on Si(OOl) and Ge(OOl), 
as well as the attractive energy of two vacancies placed in adjacent dimer rows. In the case 
of Si(OOl), a similar calculation was performed by Wang and coworkers |32| . and has served 
as useful reference for many subsequent investigations. In contrast, we have found no similar 
studies for the formation energies of defects on Ge(OOl). Here we perform such a study 
using the Tersoff potential |20| . and identify the differences between the vacancies on the two 
surfaces, Si(OOl) and Ge(OOl). We denote the formation energy of a DV complex by u, and 
that of an isolated pair of complexes lying in adjacent dimer rows (called a bound pair form 
here on) by u p . The binding energy of the pair is defined 3 as the difference between the pair's 
formation energy and the formation energies of two isolated DVs: 

w = Up — 2u. (12) 

We have calculated the formation and binding energies for different vacancy complexes, and 
have listed the results in Table 01 

We shall first discuss the formation energy u, and compare our results with those of 
Wang et al. 32 , which were obtained using density- functional calculations coupled with an 
empirical potential. The advantages and disadvantages of such a coupling scheme are clearly 
discussed in the original paper |H2j. In the case of 1-DV vacancy on Si(001), we find the 
formation energy of single DVs u ~ 0.19 eV, in good agreement with the value of 0.22 eV 
reported in Ref. |32j . For the other vacancy complexes, the agreement with Ref. (321 is poor 
because the empirical model |19j does not accurately capture the relaxations of the under- 
coordinated second-layer atoms that are present in the 2-DV and 1+2-DV complexes. Due 
to this limitation, the Tersoff potential JH] predicts that the 2-DV complex is unstable with 
respect to the formation of two isolated single DVs, «2-dv > 2ui_dv- We find however 
that the empirical potential is still able to predict that a 2-DV and a 1-DV on the same 
dimer row can bind to form a 1+2-DV complex (i.e. ui + 2dv < ^2-dv + ^i-dv)> though the 
binding energy for the 1+2DV complex is about half the value obtained from [22] . In the 
case of Ge(001) surface, the potential 2QJ fares better than the silicon parametrization |19j in 
describing the stability of DV complexes. Both the 2-DV and 1+2-DV vacancies on Ge(001) 
are stable against breaking into isolated smaller complexes, as indicated by their respective 
formation energies in Table 01 This result is in agreement with the observations of Yang et 
al., who report that 1+2-DVs and 2-DVs are predominant on the Ge(001) surface |18j . 

We now turn to a discussion of the binding energies of DV pairs w given in Table 
There are two main issues that complicate the comparison of theoretical results in TableOlwith 
the experimentally determined binding energies. First, we note that the binding energy of an 
isolated bound pair computed here and in Ref. are defined in a different manner than 
what is usually obtained from experiments. Experiments measure the binding energy when 

3 With this definition, the negative sign is retained to indicate the stability of the bound pair. 
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DV complex 


u 


(eV) 


w (eV) 


(meV/A) 


A (meV/A) 




1-DV 


0.186 


(0.22)* 


-0.281 (-0.1)1" 


-12.37 


-20.56 


Si 


2-DV 


0.505 


(0.33)* 


-0.106 


51.94 


51.57 




1+2-DV 


0.619 


(0.42)* 


-0.322 


38.66 


30.56 




1-DV 


0.549 




-0.239 


38.75 


32.48 


Ge 


2-DV 


0.699 




-0.136 


70.37 


69.61 




1+2-DV 


1.218 




-0.330 


110.00 


104.64 



* Ref.jSU 
t Ref.QI] 

Table 3: Formation energies u of vacancy complexes and binding energies w corresponding to 
a bound pair of vacancy complexes situated in adjacent dimer rows. The last two columns 
give the line formation energies estimated using formula (|13|) . as well as the ones obtained 
from total energy calculations (Table |2J). The results of Wang et al. 32 and Weakliem et al. 
|llj for the case of silicon are indicated in parentheses. 

the DVs are part of a vacancy line, while calculations address an isolated pair of vacancies. 
We have estimated that the binding of two vacancies (on an average) within a VL is up to 
~ 40% stronger than the binding within an isolated bound pair. Secondly, the experiments 
report the binding energies for DVs on Si(001) with a certain amount of Ge coverage ^Uj or 
metal contamination |13|I14|. as opposed to clean surfaces. Given these two concerns related to 
experimental conditions, as well as the use of an empirical potential in this study, discrepancies 
between the binding energies w listed in Table |3] and the ones reported experimentally are 
expected. However, we have found that there is reasonable agreement with the experimental 
data. For the single DVs on a Ge-covered Si(001) surface (1.5ML) Chen et al. estimate 
the binding energy in the range — 0.1 8eV < w < — 0.25eV, while our result for 1-DVs on pure 
Ge(001) is -0.24 eV and -0.28 eV on pure Si(001) (refer to Table©. On the Ni-contaminated 
Si(001) surfaces, Zandvliet and coworkers |14j report an average attractive interaction per unit 
length of —0.11 eV/a ; since the distance between two neighboring dimer rows is 2a, we infer 
that the binding energy of a pair within a VL is —0.22 eV. Chang et al. report a similar 
figure for the binding energy (—0.24 eV) of DVs on Ag-contaminated Si(001) surfaces |13j . 
On both Ni- and Ag- contaminated Si(001) the vacancy lines contain almost equal proportions 
of 2-DVs and 1+2-DVs ^3 El EH and we note that the binding energies reported in [T1HT3*] 
lie between the calculated (Table|3J) binding energy of a 2-DV pair (« —0.11 eV) and that of 
an 1+2-DV pair (» -0.32 eV) on Si(001). 

The vacancy formation energy u of a vacancy complex and the binding energy w of a 
pair determine the formation energy of the vacancy lines to a large extent. Chen et al. ^0] 
and Weakliem et al. jllj have found that interactions between single DVs that belong to 
different dimer rows which are not adjacent can be neglected. This approximation implies 
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that the DVs that make up a vacancy line can be treated as a 1-dimensional system with 
nearest-neighbor interactions. For such a system, the formation energy per unit length can 
be written as: 

u + w , . 

A»_, (13) 

where 2a is the distance between two DVs in a straight vacancy line. The data in Table |3] 
shows that the formation energies given by Eq. (|13|) are in good agreement with the line 
formation energies obtained from the 2 x N surface energies (Table |2J). 

The differences between the line formation energies shown in the last two columns of 
Table |3] are due to the presence of interactions between a given vacancy with higher-order 
neighbors in the same VL, which are neglected in formula (j!3j) . While such effects could 
be considered small, note from Table that they systematically lead to a lowering of the 
formation energy of vacancy lines below the estimate given by (j!3|) . As seen in the next 
section, the DVs are more stable as part of a vacancy line than as members of a bound pair, 
which is due to the presence of interactions with higher-order neighbors. 



5.3 Experimental implications 

The Si(OOl) and Ge(OOl) have been investigated intensely due to their technological and 
fundamental importance. As noted by Zandvliet |17l I43j there are still gaps in understanding 
the Ge(OOl) surface in comparison to Si(OOl), despite the similarity of the two surfaces. One 
aspect where the two surfaces show remarkable differences is that the vacancy lines appear 
easier to form on Si(OOl) rather than on Ge(OOl), as illustrated by the much larger number 
of reported observations in the case of silicon. To our knowledge, the only observation of 
vacancy lines on Ge(OOl) is for the case of Bi deposition (1ML) |45| . No vacancy lines have 
so far been observed on clean Ge(OOl) ^1 El ESI IB]) while two different groups report the 
observation of vacancy lines on the clean Si(OOl) surface |26 [ I24[ 125]. 

We look for the differences between DVs on clean Si(OOl) and Ge(OOl) by comparing 
their formation and binding energies in the presence of strain. Focusing on single DVs, we 
compare the formation energy per vacancy for isolated vacancies, bound pairs and vacancy 
lines under biaxial strain. This comparison is performed separately for the case of Si and Ge 
surfaces. The results plotted in Fig. |1] which show that the formation energy per DV, to a 
good approximation, is a linear function of strain for strains in the range -2%-2%. In contrast, 
we have found that the binding energy w = u p — 2u exhibits strong nonlinear character in 
this strain range; the strains at which w remains linear is limited to a much smaller range, 
±0.2%. Compressive strains lower the formation energy per DV both for the case of Si and 
Ge, consistent with the previous results on vacancy lines (i.e. with the positive sign of the 
biaxial strain sensitivities shown in Table |2j We can now identify the difference between DVs 
on Si(001) and Ge(001). From Table |3] note that the formation energy u of a single DV is 
positive both for Ge(001) and Si(001), but in the case of Ge(001) it is about 3 times higher 
than that of an isolated DV on Si(001). Therefore, in the absence of any kinetic effects, 
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(a) 1-DV'son Si(001) 




-0.5 - 



Figure 4: Formation energy (eV) per 1-DV vacancy for an isolated DV (u), a bound pair 
(0.5u p ), and a 1-DV vacancy line (2aA) as function of biaxial strain e. The upper and 
lower figures represent the case of Si(001) and Ge(001), respectively. Note that a and A are 
the surface lattice constant and VL formation energy, both of which depend on e. In each 
figure, the stability of DVs increases from the isolated DV to the vacancy lines, at given 
strain. Compressive strains lower the formation energy per DV both in the case of Si(001) 
and Ge(001), though in the case of Ge larger compressive strains must be applied in order to 
achieve negative formation energies. 
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a compressive strain of —0.5% would be required to create stable DVs (u < 0) on Si(001), 
which is readily achievable as result of quenching |24l 125 j . In contrast, the critical value of 
compression needed to have single DVs on Ge(001) is four times larger, —1.9%. As seen in 
Fig. |1J when DVs are part of a bound pair the compression at which the formation energy per 
DV becomes negative decreases due to the attractive interaction between DVs: for Si(001) 
the critical compression becomes —0.1%, while for DVs on Ge(001) it remains rather large, 
-1.5%. 

Based on these results (Fig. 01, we propose that organization of the DVs into vacancy 
lines on Ge(001) could be observed provided that the surface is under compressive strains 
of ~ 2%. Biaxial strains are somewhat difficult to apply and control accurately in practice. 
A simple way to provide some compression would be quenching, as done in |24| I25j . As far 
as lowering the formation energy of DVs is concerned, uniaxial strains along the dimer row 
direction are more efficient than biaxial deformations. Therefore, the use of ultraflat, single- 
domain surfaces under uniaxial compression along the dimer row can lead to observations of 
vacancies organizing into lines on Ge(001). 

6 Conclusions 

In conclusion, we have performed an investigation of the VLs and DVs on strained Si(001) 
and Ge(001). We have found that the dimer- vacancies on the two surfaces exhibit differences 
in terms of the formation energy, their strain sensitivity, and the binding energy of two DVs 
in adjacent dimer rows. With the present model of atomic interactions |19[ 120] . each of 
these energetic parameters appear to favor the self-organization of DVs into vacancy lines on 
Si(001), but not on the Ge(001) surface. 

The most important difference concerns the formation energy of single DVs, which was 
found to be 0.55 eV for Ge(001), but only 0.19 eV in the case of Si(001). The formation 
energy difference between single DVs on Si and Ge surfaces seems robust with respect to the 
empirical potential used: we have recalculated these energies with the Stillinger- Weber model 
|46[ I47| . and found similar results, though the difference between the two formation energies 
predicted by the Stillinger- Weber potential was larger. This finding cannot be explained by 
scaling the DV formation energies with the bulk cohesive energy (or melting temperature) of 
Si and Ge. Rather, it could likely be due to subtle details in the relaxation of the atoms that 
surround the vacancy, particularly to the interplay between the lateral and vertical relaxation; 
further studies are underway in order to elucidate this point. While it is well known that the 
stability of DVs on (001) surfaces is strongly enhanced by kinetic effects (i.e. two atoms are 
unlikely to find the DV and fill it at the same time [IB]); it is conceivable that the large 
gap between the formation energies of DVs on Si(001) and Ge(001) determines their different 
physical behavior. 

Based on studies of how the formation energy varies with external strains, we have 
proposed that the self-assembly of vacancies on clean Ge(001) would be possible if compressive 
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strains are applied either biaxially or along the dimer row direction. Pending an increase in the 
available computational power and methodologies, the results of this article warrant future 
studies using higher-level total-energy methods such as tight-binding or density functional 
calculations. The use of such methods is expected to bring quantitative improvements of the 
energetic properties computed here. Furthermore, higher-level methods can also address other 
important aspects of the self-assembly of vacancies such as the energetic barriers associated 
with the creation and diffusion of DVs on Si(OOl) and Ge(OOl). Future investigations regarding 
the kinetics of dimer vacancies would be invaluable for fully elucidating the differences between 
Si(OOl) and Ge(OOl). 
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